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ABSTRACT 

Frequency  domain  methods  are  siployed  for  the  design  and 
analysis  of  anti-roll  tank  systems.  A  waveslope  spectrum  of 
the  ocean  is  developed  based  upon  the  trochoidal  model  of 
ocean  wa/es  and  the  Pierson-Moskowitz  ocean  wave  power 
spectra.  The  response  of  a  typical  ship  to  the  waveslope 
spactrum  is  studied  with  the  aid  of  a  numerical  example. 

FORTRAN  subroutines  for  the  design  of  passive  and 
pump-activated  U-tube  anti-roll  tanks  are  developed,  and  the 
effects  of  a  tank  installation  on  the  ship  parameters  is 
observed.  Many  controller  design  methods  are  investigated, 
and  final  design  is  accomplished  by  exhaustive  search  of 
closed-loop  frequency  responses.  Frequency  responses  of 
passive  and  activated  tank  stabilized  ships  are  compared 
with  that  of  an  unmodified  ship  to  illustrate  the 
effectiveness  of  design.  The  effects  of  vessel  displacement 
variations  due  to  consumption  of  fuel  and  stores  on  tank 
effectiveness  and  controller  design  are  discussed. 
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I.   FORMULATION  OF  THE  PROBLEM 


A  ship  on  the  ocean  is  in  constant  motion.  Ideally,  the 
only  motions  would  be  translation  in  the  forward  direction 
and  rotation  about  the  vertical  axis,  associated  with 
propelling  and  maneuvering  the  ship.  Realistically,  there 
are  perturbation  motions  in  addition  to  the  desired  motions 
in  the  form  of  translations  along  and  rotations  about  the 
longitudinal ,  transverse,  and  vertical  axes  of  the  ship. 
Tha  purposs  of  this  thesis  was  to  study  the  rotation  about 
the  longitudinal  axis,  known  as  roll,  and  arrive  at  a  scheme 
for  controlling  it.  A  procedure  for  designing  the 
installation  of  a  U-tube  anti-roll  tank  system  was 
consolidated  from  several  sources  [Chadwick,  1950  j  [Chadwick 
and  Klottar,  1953  [Biagoveshchensky,  1 962j I  Vachanaratana, 
1973J  |_Ziegler ,  1975  and  presented  in  the  form  of  a  FORTRAN 
subroutine.  Many  controller  design  methods  were 
investigated,  and  the  final  design  was  accomplished  by 
exhaustive  search  of  closed-loop  frequency  responses. 


A.   WAVESL3PE  SPECTRUM  OF  THE  OCEAN. 


Before  a  meaningful  analysis  of  ship  motions  in  water 
could  be  conducted,  the  motions  of  the  water  itself  had  to 
be  understood.  Because  a  ship  tends  to  ride  normal  to  the 
surface  of  the  water,  the  most  significant  physical  property 
of  waves  with  respect  to  ship  motions  is  the  waveslope,  \J,  . 
At  the  crest  and  in  the  trough  of  a  wave,  the  waveslope  is 
zero,  while  at  some  point  nearly  midway  between  those 
extremes,  the  waveslope  is  maximum. 

Several  theories  of  wave  motion  have  been  advanced, 
including  the  simple  harmonic  wave;   the  Stoke's  Wave,  which 


is  a  harnonic  wave  with  a  second  harmonic  added;  and  the 
trochoidal  wave.  Geometrically,  a  trochoidal  curve  is 
constructed  by  rolling  a  circla  of  radius  R  under  a 
horizontal  line,  and  tracing  the  locus  of  a  point  a  radial 
distance  r<R  from  the  center  of  the  circle.  The  trochoidal 
curve  may  be  steep  or  shallow  depending  on  the  ratio  of  r:R, 
and  corresponds  very  closely  to  actual  observed  ocean  waves, 
where  2  rr  R  =  L  =  wavelength,  and  2r  =  h  -  trough-to-crest 
wave  height. 

Ocean  waves  are  primarily  caused  by  the  transfer  of 
kinetic  eiergy  from  the  wind,  which  is  the  result  of  the 
transfer  of  thermal  energy  in  weather  phenomena.  As  the 
wind  blows  over  a  flat  surface  of  water,  the  water  particles 
at  the  air-sea  interface  are  excited  by  frictional  forces, 
and  tend  to  move  downwind.  The  freedom  of  motion  of  these 
particles  is  limited  by  the  fluidic  properties  of  water,  and 
thus  they  tend  to  accumulate  in  successive  ridges,  or 
ripples,  on  the  sea  surface.  Once  these  ridges  are  formed, 
further  wave  buildup  continues  by  a  different  mechanism. 
On  the  upwind  side  of  a  ridge,  a  small  local  high  pressure 
region  develops  which  exerts  a  downward  force  on  the  water 
surface,  forming  the  wave  trough.  On  the  downwind  side  of  a 
ridge,  a  corresponding  small  local  low  pressure  region 
develops  which  draws  the  water  surface  upward,  forming  the 
wave  crest.  The  wave  continues  to  grow  as  long  as  the 
wind  continues  to  supply  energy  in  excess  of  that  consumed 
by  the  internal  friction  of  the  water  in  the  wave,  which  it 
does  until  the  wave  velocity  equals  the  wind  velocity. 

Three  key  factors  govern  the  development  of  wave 
systems.  They  are  the  initial  difference  between  the 
velocities  of  the  wind  and  any  existing  wave  systems;  the 
length  of  time  (duration)  that  the  wind  continues  at  a  given 
velocity;  and  the  fetch,  or  distance  of  relatively  uniform 
depth   water   over  which   the   wave   will  run   before   being 


altered  bf  a  change  of  depth  or  destroyed  by  encountering  a 
land  mass.  (The  passage  of  a  vessel  through  an  ocean  wave 
system  does  little  to  alter  the  growth  or  shape  of  the 
system. ) 

The  rate  of  increasing  wave  height  is  a  function  of  the 
initial  difference  between  the  velocities  of  the  wind  and 
the  existing  wave  system.  A  sudden  change  of  wind  velocity, 
such  as  that  associated  with  the  passage  of  a  front,  can 
cause  the  waves  to  increase  at  a  rate  of  one  or  two  feet  per 
minute.  There  is  a  physical  limit  to  the  increase  in  wave 
height,  and  it  is  defined  in  terms  of  the  wave  height  to 
wave  length  ratio,  h:L=1:7.  Near  this  ratio  limit,  excess 
wind  energy  blows  off  the  tops  of  the  waves,  forming 
turbulence  known  as  "white  caps." 

The  rate  of  increasing  wave  length  generally  lags  behind 
the  rate  of  increasing  wave  height.  Thus  in  regions  of 
short  fetch,  a  sudden  change  of  wind  velocity  may  never 
produce  waves  with  long  wavelengths  due  to  the  limiting 
ratio.  In  contrast,  on  the  open  ocean  where  there  is  a  long 
fetch,  a  change  of  wind  velocity  can  produce  long  waves  of 
relatively  low  height.  For  waves  of  length  L>100  feet,  there 
is  a  diminishing  likelihood  of  ever  approaching  the  maximum 
value  of  the  ratio.  In  practice,  for  wavelengths  between 
500  and  500  feet,  the  ratio  rarely  exceeds  1:20;  for  the 
longer  ocean  waves  (of  length  over  1000  feet) ,  the  ratio 
drops  to  1 :  50. 

The  above  discussion  of  wave  development  only  hints  at 
the  complexity  of  attempting  to  define  the  waveslope 
spectrum  of  the  ocean.  Conferences  have  been  held,  and  many 
books  have  been  written  on  the  subject  of  ocean  spectra. 
Conference,  Easton,  Maryland,  1961J  [_Inoue,  1966  J 
Comstock,  1967~]  TGillaer,  1969~|  Tziegler,  1975~]  Many  ocean 
spectra  plot  relative  wave  energy  as  a  function  of  frequency 
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with  "sea  state"  as  a  parameter.  Others  define  wind 
velocity  and  duration  as  parameters.  Still  others  define 
"significant  wave  height"  as  a  parameter,  with  many 
definitions  of  "significant  wave  height"  itself.  Spectra 
have  been  constructed  for  given  locations  citing  mean  wave 
height,  or  mean  squared  wave  height,  or  RMS  wave  height  as 
parameters . 

Two  other  factors  which  affect  the  effective  ocean 
spectra  are  vessel  speed  and  angle  of  encounter.  To 
illustrate  these  effects  on  the  ocean  spectrum,  consider  a 
ship  proceeding  with  velocity  identical  to  that  of  a 
single-frequency  wave  train.  If  the  ship  heads  directly 
into  the  approaching  wave  train,  the  apparent  frequency  of 
the  wave  doubles,  but  the  roll  coupling  with  the  ship  goes 
to  zero.  When  the  ship  travels  parallel  to  the  wave  train, 
which  is  to  say  the  ship  is  in  a  beam  sea,  the  apparent 
frequency  of  the  wave  is  the  same  as  the  actual  frequency  of 
the  wave,  and  the  roll  coupling  is  maximum.  When  the  ship 
heads  in  the  same  direction  as  the  wave  train  ("following 
seas"  coniition) ,  the  apparent  frequency  of  the  wave  drops 
to  zero,  as  does  the  roll   coupling.  If  B  denotes  the  ship's 


o 


angle  of  encounter  with  the  wave  train,   where  B=0   for  a 

0  o 

following  sea,  3=90  for  a  beam  sea,  and  B=180  for  a  head 
sea,  then  the  apparent  waveslope  is  related  to  the  actual 
waveslope,  and  the  apparent  frequency  of  the  sea  is  related 
to    the    actual    frequency    by    the   expressions   [Ziegler,    1975  |: 

<£(B)    =      +  sin    B  (1) 

max 

w  (B)    =      co  (1.-    co       u  /g)    cos    B  (2) 

sea  sea   ship 

a  =    forward   velocity   of   the   ship 

ship 

g  =   acceleration   due   to   gravity 

Since  the   ultimate      reason   for    this    study    was      to   aid   in 
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the   desigi      of   a      general-purpose   anti-roll      stabilizer,    the 

ocean      spectrum      chosen      vras    a  wave      height      power      density 

spectrum,    with      sea  state   and  significant   wave      height    (1/3 

highest   waves)      as   parameters.  The    waveslope      spectrum    was 

then   computed      on   the    basis    of  the    trochoidal    model      of    the 

ocean,    where 

r     =  1/f  =  2  jr  /  cj  seconds                   (3) 

2 

L     =  5.118  T  feet                      (4) 

\j/           =    (h/L)    x   57.3  degrees     (first    approximation) 

max  (5) 


B.       SHIP    ROLLING    IN    WAVES. 

A  ship  on  the  ocean  moves,  in  general,  with  six  degrees 
of  freedom  relative  to  the  coordinate  axis  system.  By 
convention,  the  positive  x-direction  is  forward  along  the 
longitudinal  centerline  of  the  ship,  the  positive 
y-direction  is  horizontally  outward  on  the  starboard  beam, 
and  the  positive  z-direction  is  vertically  downward.  The 
center  of  the  coordinate  axis  system  is  located  at  the 
vessel  center  of  gravity.  Rotation  about  the  x-axis  is 
known  as  roll,  with  positive  roll  angle,  <*>  ,  defined  as  an 
inclination  to  starboard.  Roll  can  be  considered  separately 
from  the  other  degrees  of  freedom,  which  reduces  the  study 
of  ship  motions  to  a  two-dimensional  problem  with  one  degree 
of   freedom. 

1 .      Definitions   and   Concepts   of    Naval    Architecture. 

The  following  definitions  and  concepts  of  naval 
architecture  jComstock,  196 7j  [Gillmer,  1969J  are  given  for 
the  benefit  of  the  raader,  and  are  illustrated  in  Figure 
1-1. 

a.       Displacement. 

The   displacement   of      a    vessel,  A ,    is      defined    as 
the    weight    of      the   volume   of    water   displaced      by    the   vessel. 
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G  =  Center  of  Gravity 

B   =  Center  of  Buoyancy 

M  =  Metacenter 

GM  =  Metacenter  Height 

GZ  =  Righting  Arm 

4>   =  Ship  Roll  Angle  Relative  to  Earth  Reference 

WL  =  Water  Line 


Figure  1-1.   Definitions  and  Concepts  of  Naval  Architecture 
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Due  to  the  Principle  of  Archimedes,  the  displacement  equals 
the  weight  of  the  vessel  itself.  Displacement  is  always 
given  in  long  tons,  where  one  long  ton  equals  2240  U.S. 
pounds.  The  weight  of  the  water  displaced  equals  the  weight 
of  the  vessel,  independent  of  the  density  of  the  water.  The 
volume  of  water  displaced,  howevsr,  is  a  function  of  the 
density  of  the  water.  For  purposes  of  most  calculations, 
the  density  of  fresh  water  is  given  as  36  cubic  feet  per 
ton,  whereas  the  density  of  salt  water  is  given  as  35  cubic 
feet  per  ton. 

b.  Center  of  Gravity. 

The  center  of  gravity,  G,  is  defined  as  the 
point  at  which  the  summation  of  the  moments  of  all  the 
weights  ia  a  vessel  with  reference  to  any  axis  passing 
through  that  point  equals  zero.  The  center  of  gravity  of 
most  warships  of  usual  form  is  generally  located  near  the 
designer's  water  line,  in  or  near  the  plane  of  the  midship 
section.  For  purposes  of  buoyancy  and  stability  analysis, 
the  entire  weight  of  a  vessel  is  considered  to  be 
concentrated  at,  and  the  force  of  the  weight  acts  vertically 
downward  through,  the  center  of  gravity. 

c.  Center   of   Buoyancy. 

The  center  of  buoyancy,  B,  is  defined  as  the 
geometric  center  of  the  underwater  form  of  the  vessel.  The 
force  of  buoyancy,  equal  in  magnitude  to  the  force  of  the 
displacemsnt  of  the  vassel,  acts  vertically  upward  through 
the  center  of  buoyancy.  When  a  vessel  is  in  static 
equilibrium,  the  center  of  buoyancy  is  located  directly 
under  the  center  of  gravity,  such  that  the  force  vectors  of 
displacement  and  buoyancy  are  equal  and  opposite.  Due  to 
this  definition  of  the  center  of  buoyancy,  as  the  condition 
of  a  vessel  varies  from  that  of  static  equilibrium,  the 
location  of  B  varies  also. 
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d.   Metacenter. 

As  a  ship  inclines,  the  shape  of  the  underwater 
form  changes,  and  the  center  of  buoyancy  migrates  to  remain 
at  the  geometric  center.  The  migrations  form  a  locus  which 
is  dependent  upon  and  directly  related  to  the  shape  of  the 
ship's  hull.  This  locus  is  generally  elliptic  and  as  such 
cannot  have  a  fixed  center,  but  for  an  infinitesimal 
increment  of  the  locus  the  instantaneous  center,  called  a 
metacenter,  can  be  located.  The  locus  of  metacenters  is 
geometrically  an  evolute,  and  is  called  the  metacentric. 
For  small  angles  of  inclination,  up  to  10  for  most  ships, 
the  metacenters  are  co-located  at  a  point  known  as  the 
Metacenter,  M. 

e.  Metacentric  Height. 

When  a  ship  is  upright,  that  is,  when  the  angle 
of  inclination  <t>  -  0  ,  the  centers  M,  G,  and  B  lie  in  a 
straight  srertical  line  over  the  keel.  The  distance  between 
G  and  M  is  known  as  the  metacentric  height,  GM.  If  M  is 
above  G,  3M  is  defined  as  being  positive. 

f.  Righting  Arm. 

When  a  ship  is  inclined,  the  center  of  buoyancy 
is  no  longer  on  the  original  vertical  centerline.  A  point  Z 
is  described  by  the  intersection  of  a  line  drawn  vertically 
upward  from  B  with  a  line  drawn  horizontally  away  from  G. 
The  distance  GZ  is  known  as  the  righting  arm,  and  for  small 
angles   of   inclination    <t>  , 

GZ    =   GM    sin  <t>        =    GM  4>  (6) 

g.  Righting  moment. 

The  righting  moment,  Mr  ,  is  the  product  of  the 
force    (displacement)    and  the   righting   arm: 

Mr    =    A  GZ    =    A  GM  <t>  (7) 

h.       Metacentric   Height   and   Stability. 

The  metacentric  height  is  a  measure  of  the 
statical  stability  of  a  vessel.  A  large  positive  GM 
corresponds    with    a   pair      of    high-f reguency   complex   conjugate 
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poles  in  the  left  half  of  the  s-plane,  far  removed  from  the 
j  u;  -axis.  A  small  positive  GM  corresponds  with  a  pair  of 
low-frequency  complex  conjugate  poles  in  the  left  half  of 
the  s-plane,  close  to  the  j«  -axis.  If  GM  =  0,  then  there 
are  two  poles  at  the  origin  of  the  s-plane.  A  negative  GM 
puts   poles    in   the  right    half   of   the   s-plane. 

2.   Equations  of  Motion. 

a.   Natural  Period  of  Roll. 

In  still  water,  a  ship  can  be  made  to  roll  by 
applying  and  removing  external  heeling  moments.  In  the 
inclined  position,  the  righting  moment  is  equal  and  opposite 
to  the  heeling  moment,  and  the  potential  energy  of  the  ship 
is  equal  to  the  work  done  by  the  external  force  causing  the 
inclination.  When  the  external  moment  is  removed,  the 
righting  moment  produces  rotation  of  the  ship  towards  the 
upright  position.  The  potential  energy  is  converted  to 
kinetic  energy  such  that  when  the  ship  is  in  the  upright 
position,  assuming  unresisted  roll,  all  the  potential  energy 
has  been  converted  to  kinetic  energy.  As  a  result,  the  ship 
will  continue  to  rotate  past  the  upright  position  until  all 
kinetic  energy  has  been  reconverted  to  potential  energy,  at 
the  opposite  limit  of  roll.  Assuming  no  loss  due  to 
friction,  the  ship  would  continue  to  oscillate,  or  roll, 
indefinitely  with  constant  amplitude.  The  period  of  these 
oscillations  can  be  found  from  simple  harmonic  motion 
relationships  to  be  I  Comstock,  1967 J 
T   =  (27rk)/N/iGF 

=  (1 .108k)  /v/gF 


where 


~   (0.44b)  />/ GM  (8) 


2 

g  -   acceleration  due  to  gravity,  ft/sec 

k  =  radius  of  gyration  of  mass  of  ship,  feet 
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b  -   extreme  beam  of  ship,  feet 

b.  Roll  by  Wave  Action. 

The  cause  of  ship  rolling  in  a  seaway  is 
primarily  the  unbalanced  moments  which  result  from  the 
center  of  buoyancy  shifting  in  response  to  the  change  of  the 
waterplane,  and  thus  the  shape  of  the  underwater  portion  of 
th9  hull,  as  waves  pass  under  the  ship.  A  ship  will  usually 
attempt  to  ride  normal  to  the  surface  of  a  wave.  A  ship 
also  moves  in  a  circular  orbit  as  do  the  particles  of  water 
in  a  wave.  The  orbits  induce  centrifugal  forces,  which  are 
opposed  by  the  dynamic  forces  in  water. 

Rolling  in  waves  is  a  forced  oscillation,  and  the  ship 
will  tend  to  oscillate  at  or  near  the  frequency  of  the 
waves.  If  the  waves  are  not  uniform,  which  is  usually  the 
case,  a  ship  will  attempt  to  oscillate  with  its  natural 
frequency.  The  overall  roll,  therefore,  is  a  combination  of 
tha  ship»s  natural  and  the  wave's  periods,  with  the  wave 
influence  being  generally  more  predominant. 

c.  Equations  of  Motion. 

A  ship  can  be  modeled  as  a  second-order, 
spring-mass-damper  system.  [Chadwick,  1950~]  Fcomstocic,  1967"] 
The  spring  constant,  or  static  righting  coefficient,  was 
discussed  earlier  in  this  paper  in  the  definition  of  the 
righting  moment.      The   static   righting   coefficient  is 

K      =AGM  (9) 

s 

The   inertia      term,    or   roll     inertia   coefficient,      is    the 

product  of  the      mass   of   the   ship   and   the      radius   of   gyration 

of   the    mass    of  the   ship   squared: 

2 
J      =    (  A/g)    k  (10) 

s 

The  ridius   of  gyration   is  found   by  timing   the  natual 
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period  of  roll  in  protected  (still)  water,  and  massaging  the 
equation  for  roll  period  as  follows: 
k   =  <T+>/gG¥  )/2* 

=  (T^^/sT  )/1.108  (11) 

The      roll      damping      coefficient  is     derived      from      the 

observation      that   the      damping   ratio  of  a      typical   ship      is 
approximately 

f  =   0.1/tt  =    0.0318  (12) 

ship 

By      working   backwards      through      the   second-order     system 

relationships,    the   roll   damping   coefficient   was    found   to    be 

2 
B         =      0.1     (I.   A  Gfl    )     /  7T  (13) 

s  ♦ 

On  the  forcing  function  side  of  the  equation  of  motion, 
following  the  classical  usage  of  Froude  and  Krylov, 
[chadwick  and  Klotter,  1953  |  the  ship-sea  coupling 
coefficient  was  considered  to  be  approximately  equal  to  the 
static  righting  coefficient.  Therefore,  the  equation  of 
motion  for  a  ship  in  a  seaway  is: 

J  £  +  B  <*>   ♦  K  <*>   =   K    \p  (14) 

s      s       s       s 

Dividing   through  by   the  roll   inertia  coefficient   and 

performing   the   LaPlace  Transform   operation,   the   system 

transfer  function  was  obtained: 

fs      ♦     (B  /J  )s   +    (K    /J   )")  $(s)     =    (K   /J    )     *(s)  (15) 

*-  ss  ss-J  ss 

2  2  2 

\~S        +      (2fU)S     ♦        ">  *(S)      »     U      *(S)  <16> 

L-  S     S  S  J  S 

2 

ct>   (S)  o} 

_S_ 


^2  2 

^  (s)  s      +     (2  f  ci>  )s   +   u 

s   s  s 
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(K  /J  ) 

s  s (17) 

2 
s   ♦  (B  /J  )  S  ♦  (K  /J  ) 

s   s        s   «. 

It  was  useful  in  the  analysis  to  normalize  the  equations 
such  that  the  ship's  natural  frequency  was  unity.  This  was 
accomplished  in  the  following  manner: 

Starting  with  the  general  equation  of  motion,  divide 
through  by  the  static  righting  coefficient: 

(J    /K    )   <t>  +     (B    /K    )  <t>  +  <p    =    <£  (18) 

s     s  s     s 


Let    t'     =    to  t    =  f/K    /J    ) t 
s  s      s 

Define    d  <t>  /dt  =    <t> 

d<f>/dt'  =    <t>' 

2  2 

d  <t>  /dt  ■  =    <*>" 


Then      <t>  -    UK  /J      )  <t>' 
s     s 

$  =    (K   /J   )    </>" 
s      s 

The   equation   was   then    rewritten   as 

<*>"♦     (B    /   K   J      )   <*>'  ♦   <t>    -    *  (19) 

s        s    s 

The  LaPlace  Transform  operation  was  performed,   and  the 

system  transfer  function  (normalized)  was  obtained: 

Ts   ♦  (B  /  K  J   )s  +  1~]  *  (s)  *  y(s)  (20) 

L-  S     S  S  J 

*  (s)  1 


2  (21) 


*  (s)  s   ♦  (B  /  K  J  )  s  +  1 

s   s  s 


3.   Frequency  Response. 


In  order  to  study  the  frequency  response  of  a  vessel 
to  the  waveslope  spectrum  of  the  ocean,  it  was  necessary  to 
introduce  a   numerical  example.    The  example   chosen  was   a 
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well  used  example  [Blagoveshchensky,  1962 j,  which  was: 
L  =   49.1  meters      (161.0  feet) 
b  ■   S.O  meters       (  29.5  feet) 
A  =  936.0  metric  tons  (921.2  long  tons) 
GM  =  0.73  meter       (   2.4  feet) 
T.  =   8.63  seconds 
The  radius  of  gyration  of  the   mass  of  the  ship,  and  the 
coefficients   of  the   equation  of   motion  were  computed   as 
follows : 


k  =  (8.637(9781)  (0.73)  )  /2*  =  3.68  meters 

2  2 

J    =    (936/9. 81) (3.68)       =    1289.021    tons-meters-sec 

s 

2 
B    =    (0.  1)  (8.63)  (936)  (0.73)  /  *     =59.7    tons-meters-sec 

3 

K  =  (936)  (0.73)  =  683.28  tons-meters 
s 

These   values  were   then  substituted  into  the  transfer 

function,   the  roots   of   the  characteristic  equation   were 

found,  and  the  frequency  response  was  determined. 

<t>  (s)   =         0.53 

—      ^ (22) 

vj,  (s)       s   ♦  0.046s  ♦  0.53 

roots  (poles)  at   s  =  -0.02  ±   j  0.73 

natural  frequency   u»  =  0.73 

s 

damping  ratio      f  =  0.032 

s 

The   frequency  response   is  graphically  illustrated   in 

Figure   1-2.   The  Digital   Simulation   Language   (DSL/360) 

program  which   generated  the  plot   is  listed  in   Appendix  A. 

The   zero  decibel   low  frequency   response  demonstrates   the 

natural  tendency  of  a  ship   to  ride  normal  to  the  surface  of 

a  wave.   The  resonant  peak  which  occurs  at  0.724  rad/sec  has 

a   gain  of   23.9  db,  or  a   magnification   factor  of   15.6. 

According  to  the  waveslope  spectrum   of  the  ocean,  tabulated 


20 


FREQUENCY  RESPONSE,  SHIP  WITHOUT  TRNKS 
oCLGSED  LOOP  RESPONSE,  MflG  (DB)  VS.  LOG  W 


-1.00 


-0.80 


0.80 


Figure  1-2.   Frequency  Response  of  Typica  Ship 
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in  Appendix  E,  for  a  sea  of  15  foot  significant  wave  height 
(sea  state  6)  the  maximum  waveslope  angle  at  0.724  rad/sec 
would  be  2.23  degrees.  Therefore,  the  ship  would  experience 
a  maximum  roll  angle  of  34.8  degrees,  relative  to  the 
horizon.  This  severity  of  roll  causes  crew  discomfort,  and 
limits  the  ability  of  the  ship  to  conduct  operations,  such 
as  gunfire  support,  flight  operations,  or  small  boat 
operations. 
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II.       SHIP    STABILIZATION    BY    PASSIVE    ANTI-ROLL   TANKS 

A.       DEVELOPMENT    OF    THE    EQUATIONS    OF    MOTION 

In  order  to  apply  servo  theory  to  the  problem  of  ship 
stabilization,  the  physical  system  was.  assumed  to  be  linear, 
with  lumped  parameters  and  constant  coefficients.  It  was 
also  assumed  that  the  rolling  degree  of  freedom  could  be 
considered  independently  of  the  other  degrees  of  freedom, 
such  that  the  ship  oscillates  as  if  about  a  fixed  center  of 
rotation.  The  problem  thus  became  a  two-body  problem,  the 
ship  and  the  tank-fluid,  each  with  one  degree  of  freedom. 
The  tank-fluid  was  idealized  as  a  one-dimensional  filament, 
with  provision  for  varying  cross-section  but  with  constant 
velocity  across  each  cross-section.  These  assumptions 
required  small  angles  of  ship  roll  but  allowed  for  moderate 
angles  of  tank-fluid  level,  both  being  reasonable  in  a 
stabilized    ship. 

1 .      Double   Pendulum    Analogy. 

As  discussed  earlier  in  this  paper,  a  ship  at  sea 
behaves  like  a  damped  oscillator,  and  was  modeled  as  a 
second-order  spring-mass-damper  system.  Another  appropriate 
model  which  was  considered  was  the  pendulum  model  Tchadwick, 
1950J  Tchadwick  and  Klotter,  1953 |,  with  a  spring  constant, 
an    inertia   constant,    and   a    damping   constant. 

The  stabilizer  chosen  was  a  U-tube  anti-rolling  tank 
system.  The  water  or  other  fluid  in  the  U-tube  acts  as  a 
hydraulic  pendulum.  Since  the  tank  was  to  be  contained 
within  and  wholly  supported  by  the  ship,  the  ship-tank 
system  thit  resulted  was  a  pendulum  supporting  a  pendulum, 
or    a   double    pendulum.       In    general,    this   kind   of    system  leads 
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to      two   second-order      differential      equations    with     coupling 
terns . 

2.   Nomenclature. 

There  are  two  types  of  nomenclature  needed  here. 
These  are  the  parameters  and  variables  of  the  differential 
equations,  and  the  dummy  variables  and  parameters  which 
characterize  the  tank  geometry. 

a.   Differential  Equations. 

The  differential  equations  are  second-order  in 
nature,  and  have  forcing  functions.  The  variables  and 
parameters  are  as  follows: 

Dependent  variables  of  the  forcing  side: 

^  =  instantaneous  effective  waveslope 
Parameters  of  the  forcing  variables: 
K  =  ship-sea  coupling  coefficient 

S3 


Torques  produced  on  the  forcing  side: 

K   ^  -   disturbance  torque  due  to  waves 
ss 

Dependent  variables  of  the  homogeneous  side: 

<t>   -    ship  roll  angle  relative  to  the  horizon 

0  =  tank-fluid  level  angle  relative  to  the  ship 

Parameters  of  the  homogeneous  variables: 
Ship  self-terms: 
j»  =  roll  inertia  coef.,  unmodified  ship  (no  tanks) 

3 

B*  =  damping  coef.  of  unmodified  ship 

3 

K1  =  static  righting  coef.  of  unmodified  ship 

3 

J      =    roll   inertia   coef.,    modified  ship    (with   tanks) 

3 

3  =  damping  coef.  of  modified  ship 

3 

K  =  static  righting  coef.  of  modified  ship 
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Tank  self-terms: 

J  =  equivalent  inertia  coef.  of  tank-fluid 

t 

B'    =    skin-friction    damping   coef.    of   tank-fluid 

B      =    overall   damping   coef.    of   tank-fluid 

t 

K      =    equivalent    static  coef.    of    tank-fluid 

t 

Hitual  coupling  terms: 

J   =  mutual-coupling  inertia  coefficient 

3  t 

B   =  mutual-coupling  damping  coefficient 

3  t 

K     =    mutual-coupling   static    coefficient 

3 1 

It  is  a  feature  of  the  double  pendulum  that 

K   =  K 

St     t 

B   =  0 
-  St 

As  discussed  earlier  in  this  paper 

K   =  K 

35     S 

b.   Tank  Geometry. 

The   notation   used   in    the   derivation   which 

depends  upon   tank  geometry  is  shown   in  Figure  2-1 ,   and  is 

given  belDw. 

s    =   a  generalized   coordinate,   the   curvilinear 

distance  along  the  fluid   trajectory.   s=0  at 

the  tank-fluid  surface  in  the  starboard  tank, 

and  s=S  at  the  tank-fluid  surface  in  the  port 

o 
tank,  when  0  =0. 

R  =  lever  arm  of  tank,  the  horizontal  distance 
from  the  vertical  centerline  of  either  side 
tank  to  the  vertical  centerline  of  the  ship. 

r(s)  =  the  radial  distance  from  the   virtual  center 
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Figure  2.1.   The  Dependent  Variables  and  Nomenclature  of  Tank 
Geometry 
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of  rotation   to  any    point   on  the   fluid 

trajectory. 
P         -   density  of  fluid  in  tanks. 
D(s)  =  the  perpendicular  distance  from   the  virtual 

center  of  rotation  to  tangent  at  any  point  on 

trajectory. 
a(s)  =  angle  between  r  (s)  and  D  (s)  . 
7(s)  =  angle  between  r(s)  and  deck  of  ship. 
K  =  tank  cross-section  at  s=0  and  s=S. 

a(s)  =  tank  cross-section  at  any  point  s. 
g    =  acceleration  of  gravity. 

3.   Euler-LaGrange  Equations  of  Motion. 

Tie  two  degrse  of  freedom  system  of  ship  and 
tank-fluid  was  best  treated  by  means  of  the  Euler-LaGrange 
equations  in  generalized  coordinates.  The  generalized 
equation  is  Ichadwick  and  Klotter,  1953  | 

a  D 

♦  3  q.  =  Q  (23) 

l     l 


d 

dT 

dT 

au 

dt 

1- 

♦      dq       ♦ 
i 

1 

whsre  Q.  is  the  so-called  generalized  non-conservative  force 
in  the  ith  equation,  and  i=1,2.  The  generalized  force  is 
equal  to  the  rate  of  doing  virtual  work  against  the  external 
non-conservative  forces  when  q.  is  given  a  virtual 
displacemaat .  The  parameters  of  the  Euler-LaGrange  equation 
are: 


T   =    2  mv    /2 

U   =     2  ff'&x 

^  2 

D    =     2  3  v    /2 


=    kinetic    energy 


(24) 


=  conservative  potential  energy 
(due  to  gravity)  (25) 

=  dissipative  energy  (26) 

(due  to  skin  friction,  etc.) 


f  =  conservative  force 
/3  =  damping  parameter 
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4.   Calculation  of  the  Energy  Functions. 

Foe  each  type  of  energy  (kinetic,  conservative 
potential,  and  dissipative) ,  the  contribution  of  the  ship 
alone  and  the  contribution  of  the  tank-fluid  alone  was 
calculated,  then  the  contributions  were  added  to  arrive  at 
the  total  energy.  In  the  regions  where  the  fluid  surfaces 
travel,  namely  the  vertical  legs  of  the  U-tube,  it  was 
assumed  that  the  tanks  were  straight,  parallel,  and  of 
uniform  cross-section.  This  meant  that  a (s)  and  D(s)  were 
constants  in  these  regions,  thus  for  any  expressions  which 
depended  only  on  the  variables  a(s)  and  D(s)  ,  the  limits  of 
integration  were  fixed  at  zero  and  S,  without  approximation. 

Calculation  of  kinetic  energy: 


fl 


S  +  R  tan0 


uid  =  /***■»{  [' 


R  A 
aTsT 


i  r  (sT 

2 
♦    2 

d  R     & 

.  a{s)J 

4>D(s) 


R  tanfl 

•  2 
J»  <*>  /2 

ship       s 

.       =    T       +  T 
total      fluid    ship 


Calculation  of  potential  energy: 


U 


fluid 


3  hip 


otal 


P   g  A  R 


s)  sin>  ds 


K'  <t>  /2 
s 

D       +0 
fluid     ship 


1  - 

cos  <t> 

2 

tan  $ 

( 

cos  <t> 

30) 

(31) 

( 

32) 
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Calculation  of  dissipative  energy: 


fluid 


ship 


.2 
B*  9    /2 
t 


3'<t>   /2 

s 


D  =       D  ♦    D 

total  fluid  ship 


5.      The    Equations   of   Motion. 


(33) 

(34) 
(35) 


The    equations    of   motion   of      the    ship   and   tank    system 

come      from      equation      (23)  ,      with      q,       -    <t>  ,      and      q„      =    9    . 

Substituting    in   equations    (27)    through    (35)    and   simplifying, 

and   letting    Q      =    K    ^  and  Q„    =   0,    the    equations   of   motion    are 

»  L 

Ship:       {J    '4  *B    i  *K   *  )    ♦    (J       9  +  K    9   )    =    K   ^  (36) 

s  s  s  St  t  s 


Tank:       (J      4>  +K    <t>  )     ♦     (J    9   +B    9    *K    9    )    =    0 

st  t  t  t  t 


(37) 


Obtaining  the  system  transfer  function  from  the 
equations  of  motion  was  somewhat  more  complicated  than  it 
was  for  the  ship  without  tanks.  First,  the  equations  were 
divided  through  by  the  respective  inertia  coeficients,  and 
the    LaPlace    Transform   operation   was    performed: 

"?  o 

Ship:    fs    +(B    /J   )  s+  (K   /J    fU  (s)  ♦  [~(J      /J    )s    ♦  (K   /J    fle(s) 

u  ss  ss-1  L-     st      s  ts-1 


=     (K   /J    )    *(s) 
s     s 


(38) 


7  7 

Tank:     [~(J     /J   )s   +  (K   /J    )~U(s)*Ts    *  (B  /J    )  s+ (K  /J   )le(s) 
Lstt  ttJ  L  tt  t     t  -1 

(39) 


=    0 


4>  (s)    as      isolated    by      applying    Cramer's      Rule,    and      the 
transfer   function  was   obtained: 
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A    =    s    ♦ 


B    J    *J    B 

St        St 

r 

j  j  -j 

S     t        St 


3 

s   + 


K    J   +J    K    +B    B   -2K    J 

s   t      st      st        t    S± 

J    J    -J 

s   t      st 


B    K    +K    B 
s    t      s    t 

T~ 

J  J   -J 
s   t      st 


s    + 


K    K    -K 
s    t      t 
2 

J    J    -J 
s   t      st 


(40) 


<t>   (3)  (K    /J    )       s    ♦  (B   /J    )s+(K    /J    ) 

=  s       s t       t t       t 

*    (3)  £ 

6.      Definitions   of    Parameters. 


(41) 


The  parameters  of  the  equations  of  motion  are  given 
in  terms  of  previously  define!  parameters  and  dummy 
variables   as    follows: 


Ship    self-terms : 


2    /*  2 

J      =    J1    ♦    PAR      /  [a  (s)/A~|  fr  (s) /R~|      ds 


(42) 


B      ■    B' 
s  s 


K      =    K»    +  pgAR    /Ta  (s)/A~|  [r  (s) /r1  sin  7     ds 


(43) 


(44) 


Tank    self-terms: 


J      =    P   A    R       /     A/a  (s)       ds 

0 

B      =       B» 
t  t 


K      =      2  P  g    A    R 

t 


(45) 
(46) 
(47) 
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Mutual-coupling    terms: 

S 

J      =    p    A    R       /      D(s)  /R      as  (48) 

5t  J 

0 
B      =    0  (49) 

3  t 

K      =    K  (50) 

st         t 


B.   ANALYSIS  OF  THE  EQUATIONS  OF  MOTION. 

Having  arrived  at  a  set  of  equations  of  motion,  it  was 
necessary  to  conduct  a  brief  analysis  of  the  equations  to 
give  them  physical  significance. 

1 .  The  Self  Terms. 

Looking  first  at  the  self-terms  of  the  homogeneous 
equations,  it  was  observed  that  the  ship  acts  as  a  simple 
pendulum  if  the  tank-fluid  does  not  move,  and  that  the 
tank-fluid  acts  as  a  simple  pendulum  if  the  ship  is 
prevented  from  moving.  Each  of  these  pendulums  has  its  own 
natural  frequency.  These  two  frequencies  are  important 
factors  in  the  stabilization  problem. 

2.  The    Mutual   Terms. 

There  is  a  third  critical  frequency  associated  with 
the  ship-tank  system,  called  the  secondary  resonance 
frequency,  which  is  due  to  the  mutual-coupling  terms.  The 
mathematical  significance  of  the  effect  became  apparent  when 
the  ship  equation,  equation  (36) ,  was  rewritten  with  the 
mutual   teems    on   the  right   side  as   follows: 


J    <p  +B    <t>  +K    <J>   =    K    +   -  (J      6    +K    0    )  (51) 

3  S  S  S  St  t 
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This  aquation  shows  the  ship  as  a  simple  pendulum,  acted 
upon  by  a  net  torque,  which  is  the  torque  due  to  waves  minus 
the  torque  due  to  tank-fluid  motion.  If  the  tank-fluid  were 
made  to  oscillate  at  the  secondary  resonance  frequency,  the 
tank-fluid  torque  would  be  equal  to  the  wave  torque,  and 
thus   the  ship    would   be    influenced    by    zero  net   torque. 

C.       ANTI-ROLL    TANK    DESIGN    CRITERIA. 

The  easiest  tank  geometry  to  dasign  and  build  is  a  right 
rectangular  U-tube,  with  the  vertical  legs  having  constant 
cross-sectional  area,  and  the  horizontal  cross-leg  having 
constant  but  smaller  cross-sectional  area.  Given  this  tank 
geometry, 

a(vertical   leg)         =      A  (52) 

a(horizontal  leg)  =   a 

2  2 

=  (A  R  w  )/(g  -  h  u    )  (53) 

s  s 

D(vertical  leg)    =   R  (54) 

D(horizontal  leg)  ■   D  (55) 

The  length  of  the  trajectory  of  the  fluid  in  the  tank, 
S,  is  dependent  upon  the  lever  arm  of  the  tank,  R,  and  the 
height  of  the  fluid  column  in  a  vertical  leg  of  the  U-tube 
when  8=0°,    h, 

S  =  2  (R  +  h)  (56) 

For  maximum  effect,  the  tank  system  must  span  the 
extreme  beam  of  the  ship.  Therefore  the  lever  arm  of  the 
tank  is  half  of  the  beam,  b,  diminished  by  the  width  of  a 
vertical  lag,  w, 

R  =  (b-w)/2  (57) 

The  relationships  derived  froi  the  equations  of  motion 
which  defiae  the  tank  dimensions  and  parameters  are: 
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K    /J      =    X  (58) 

t       s 

2 
K.  /J      =     w 

*      *  t 

2    q 


/     A/a(s)     ds  (59) 

0 

J      /J      =    X 

st      t 

s 

X        =       |(D    a)/(R    A)       ds 
0 
=    2    h    ♦    2(D    a/A)  (60) 

The  first  ratio,  X ,  reflects  the  free  surface  effect  of 
the  tank-fluid,  and  should  be  a  small  number.  An  acceptable 
value    was   found  to   be   0.05    Tvachanaratana,    1973 1. 

To  be  effective  as  a  stabilizer,  the  anti-roll  tank 
system  must  be  "tuned"  to  the  natural  frequency  of  the  ship 
in   which  it    is   installed.      Therefore, 

2  2 

t  S 

K    /J       =      K    /J  (61) 

t      t  s      s 

The  third  ratio,  X,  is  dependent  upon  the  placement  of 
the  tanks  in  the  ship  relative  to  the  center  of  rotation. 
This  ratio  can  be  forced  to  zero  (J  ,  forced  to  zero)  by 
placing  the  horizontal  leg  of  the  U-tube  above  the  center  of 
rotation  by    the  amount 

-0    =    h    A/a  (62) 

There  are  no  hard  and  fast  rules  for  determining  a  value 
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of    h,    except    that    the    tanks    must      fit    within   the    hull    of    the 

ship,    and      the   fluid    in    the    tanks    should   not    be    permitted   to 

either  contact  the   top    or   completely   drain   out   of  the    bottom 

of      a   vertical     leg      of    the      tank      at      the   maximum      expected 

taak-fluid   level   angle,    0ma](  . 

Since   J        depends   upon    tank  geometry,      the   approximation 
s 

must      be  made      that   J        =   J'       ,    for      purposes   of      initiating 

s  s 

design.      Then 

K      =         X   J1 

t  s 

2 
=      2  p  g   A   8 

2 
=     P  g    1    w    (b-w)     /2  (63) 

A   reasonable    value    of   tank   width    should   be   chosen.      Then 

2 
1        =        2K/pgw     (b-w)  (64) 

Making  the  approximation  that  the  natural  frequency  of 
the  modified  ship  is  the  same  as  the  natural  frequency  of 
tha  ship  without  tanks,  and  picking  a  value  for  h,  equation 
(53)  is  solved  for  a,  and  then  equation  (62)  is  tested  for 
feasibility.  If  the  value  of  D  so  computed  cannot  be 
physically  realized,  then  a  realizable  value  of  D  must  be 
chosen.  This  merely  means  that  J  .  will  not  be  identically 
zero,    and   the   overall   analysis   will   take  more    care. 

The  tank  damping  parameter  is  derived  by  working 
backwards  through  the  second-ordar  relationships  for  the 
hydraulic  pendulum,  noting  that  the  damping  ratio  for  the 
taak-fluid  is  approximately  equal  to  the  damping  ratio  of 
tha    (unmodified)     ship: 

r  =       0.1/tt  (65) 


B         =      0.2  JK    J      / 

t  v    t    t 


(66) 


Once  the   tank  system   parameters  have   been  determined, 
the  modified  ship  parameters  are   found.   Over  the  range  s=0 
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2  2  2 

to   s=h,     r      -   R  ♦       (D  +  h-s)       ,    and      r    sin  1-   D+h-s;      over    the 

7  2  2 

range   s=h   to    s=h+R,    r        =    (R+h-s)       +    D.      ,    and   r   sin  7  =   D. 

2 
J         =      J*    +    2  P  A    R       (I    +    I    )  (67) 

s  s  12 

I         =        f(3R    +  2D    )h    ♦    2Dh      ♦    h  J/    [3R   1  (68) 

2  2 

I         =      a     (R      +3D      )/(    3AR    )  (69) 

K         =      K'    +2pg   4    R    (IM    ]  (70) 


3 


s  3         4 


2 
I        =      -  (2D    +    h    ) /(    2R   )  (71) 


T 
J. 

4 


=      Da/A  (72) 


The  design  steps  should  then  be  retraced,  removing  the 
approximations.  Normally,  after  two  or  three  iterations  of 
the  design  procedure,  the  calculated  values  converge  on  the 
actual  values,  and  a  final  design  has  been  achieved.  The 
FORTRAN  subroutine  used  in  this  iterative  design  procedure 
is   included    in  Appendix    B   to   this   thesis. 

D.       FREQUENCY    RESPONSE    OF    PASSIVE    SHIP-TANK    SYSTEM. 

The      same      numerical      example        used      in      the      frequency 

response  study     of   a   vessel      on  tha      ocean   was    used      for   the 

design    example  and      frequency    response    study    of      a   ship   with 

passive   anti-roll     tanks   to      the    waveslope     spectrum   of      the 

ocean: 

2 
J1       =      1289.021    metric  tons-meters-sec 

5 

B'       =  59.746    metric   tons-meters-sec 

5 

K»       =        683.280    metric  tons-meters 

5 

The  following  arbitrary  choices  were  made: 
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w 

5 

1 . 0      meter 

h 

= 

1 . 3      meter 

X 

S 

0.05 

Then 


p  g  =   1.025  metric  tons/meter   (salt  water  used) 

=   (0.05)  (1289.021) 

=   64.451   metric  tons-meter 

2 
(2)  (64.451)  /  (1.025)  (1)  (9-1) 

=   1.965   meter 

=   v  1 

2 
=   1.965   meter 


(1.965 )(4)  (683. 28/1289.021) 
9.81    -     (1.3) {683. 28/1239. 021 


=      0.457    meter 

Testing  equation  (62)  for  feasibility, 
-D   =    (1.3)  (1  .965)  /   0.457 
=  5.59   meter 
Since  this  tank  placement  would  probably  not  be  possible 
in  the  ship  being  studied,  a  more  realistic   value,  still 
above  the  center  of  rotation,  was  arbitrarily  chosen: 

-D   =   1.0   meter 

Then  X    is   non-zero,    as    indicated. 
J      /J      =      X 

3t         t 

=      2    \j.3    -     (0.457)  (1.965)] 
=      2.135 

To  summarize,  the  tank  dimensions  found  thus  far  were 
w   =   1.0     meter 
1   =   1.965   meter 
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2 
A   =   1.965   meter 

2 

a   =   0.457   meter 

R   =   4.0    meter 

D   =  -1.0     meter 

Given   these   tank   dimensions,    the   modified   ship 
parameters  became: 

J    =   1289.021  ♦  11.155 

2 

=   1300.176   metric  tons-meters-sec 

K    =    683.280  -   3. 121 

=        680.159      metric   tons-meter 

After  four  iterations      of   the      design      procedure,      the 

following  final  values    were   obtained: 

w        =  1.0  meter 

1        =  1.982      meter 

2 
A        =  1.982      meter 

2 

a        =  0.454      meter 

R        =  4.0  meter 

D        =  1.0  meter 

2 
J        =    1300.236      metric   tons-meters-sec 

5 

B    =    59.746   metric  tons-meters-sec 

3 

K        =      680.185      metric   tons-meters 

2 


3 


J    =   124.276   metric  tons-meters-sec 

t 

B   =    5.722   metric  tons-meters-sec 

t 

K   =   65.012   metric  tons-meters 
t 

2 
J    =   266.152   metric  tons-meters-sec 

st 

These   values  were   then  substituted   into  the   transfer 
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function,  the  roots  were  located,  and  the  frequency  response 
was  obtained. 

2 

<P   (s)  0.523    t    s      +    0.0460    s   +    0.523    ) 

"4"  3  2 


*  (s)  s      ♦    0.164    s      +    1.485    s      ♦    0.0857   s    +    0.441 

(73) 
zeroes   at        s      =      -0.0230    +    j    0.723 
pDles      at        s      =      -0.0135    ±    j    0.642 
and      s      =      -0.0684    ±    j    1.031 

The  frequency  response  is  graphically  illustrated  in 
Figure  2-2.  The  DSL/360  program  which  generated  the  plot  is 
listed  in  Appendix  3.  The  negative  decibel  low  frequency 
response  demonstrates  the  ability  of  the  passive  tank  system 
to  overcome  the  natural  tendency  of  the  ship  to  ride  normal 
to  the  waveslope,  making  it  ride  upright  with  respect  to 
earth   reference    instead. 

Comparing  this  response  with  the  response  of  the 
unmodified  ship,  given  in  Figure  1-2,  it  was  observed  that 
the  effect  of  adding  passive  anti-roll  tanks  was  to  place  a 
pair  of  complex  zeroes  at  the  ship's  own  natural  frequency 
to  reduce  the  resonant  peak,  from  24  db  at  0.7  rad/sec,  to 
two  lower  peaks  of  14  db  at  0.6  rad/sec  and  10  db  at  1.0 
rad/sec,  and  a  notch  of  -11  db  at  0.7  rad/sec.  In  sea  state 
6,  this  translates  to  maximum  ship  roll  angles  of  8.5 
degrees  at  0.6  rad/sec,  and  15.3  degrees  at  1.0  rad/sec,  a 
considerable  improvement  over  the  34.8  degree  maximum  roll 
angle   of  the    unmodified   ship. 
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FREQUENCY  RESPONSE,  PASSIVE  TANKS 

LOOP  RESPONSE,  MAG  (DB)  VS.  LOG  N 


gCLOSED 

o 


-1.00 


-0.80 


-0.60 


-0.20 


0.80 


Figure  2.2.   Frequency  Response  of  Ship  with  Passive  Tanks 
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III.   SHIP  STABILIZATION  BY  ACTIVATED  ANTI-ROLL  TANKS 

A.   MODIFIED  EQUATIONS  OF  MOTION. 

Activiting  an  anti-roll  tank  system,  mathematically, 
only  involves  adding  a  forcing  function  term  to  the  tank 
equation  of  motion.  Physically,  however,  it  can  be 
accomplished  in  basically  two  different  ways.  These  methods 
are  pneumatic  and  hydraulic.  A  pnsumatic  effector  operates 
by  alternately  creating  an  overpressure  or  a  partial  vacuum 
condition  in  the  airspace  above  the  tank-fluid  surface  in 
the  vertical  legs  of  the  U-tube,  while  an  hydraulic  effector 
generally  consists  of  an  arrangement  of  pumps  in  the 
horizontal  cross  duct.  The  effector  type  chosen  for  this 
analysis  was  a  constant  speed  controllable  pitch  propellor 
pump,  so  chosen  for  its  linear  blade  pitch  angle  vs.  flow 
characteristic  [Blagoveshchensky,  1962J. 

1 .   Nomenclature. 

In  addition  to  the  nomenclature   of  the  passive  tank 

system,  the  following  are  the  parameters  and  variables  of  an 

activated  tank  system: 

Dependant  variables  of  the  forcing  side; 

a  ~    instantaneous  pump-blade  pitch  angle  (commanded) 

Parameters  of  the  forcing  variable: 

B  =  velocity  head-loss  coefficient  of  the  pump 
P 

K  =  static  head  coefficient  of  the  pump 
p 

Torques  produced  on  the  forcing  side: 

K  a  -      B  B     -    Counter-rolling  torque  due  to  pump 
P      P 
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2.      Modification  of  Parameters. 

As    indicated      by   the   expression      for   counter-rolling 

torque,    tha      addition   of   a    pump     to    the   tank   system      has    the 

effect   of   increasing   the  overall      damping   coefficient   of    the 

tank-fluid.      All   other      parameters   of      the   homogeneous      side 

remain    unchanged.      The    pump    should      be    chosen    such      that    the 

composite   tank  damping   factor   is   0.45    Ivachanaratana,    1973  j: 

B  =    B*    ♦    B  (74) 

t  t  p 

B    /J      =    0.45  (75) 

B  =    0.45    J      -    B'  (76) 

p  t  t 

3.      Equations   of    Motion. 

Tha    equations   of   motion     for   the   activated   anti-roll 
tank   system    come    from      equation    (23)  ,    with   Q      =    K.'A    ,    Q      =    K.a 
,    q      =    <}>  ,    and  q?    =    6     I  Chadwick    and    Klotter,     1953J. 

Ship:       (J    '<$  *    B    i  f    K    <t>  )    +    (J      V  ♦    K     d  )     =    K    ^  (77) 

s  S  S  St  t  s 

Tank:       (J      £  ♦    K    <*>   )     ♦     (J    0*  ♦    B    0    +    K    0  )     =    K    a  (78) 

st  t  t  t  t  p 

Following   the      procedure    used   for   the      passive   ship-tank 
system,   <t>(s)     was    found   to    be 

$(s)       =    +    (s)     -    4>    (s)  (79) 

<t>     (s)     =    (K   /J    )fs      ♦    (B   /J    )  s    ♦     (K   /J    )"|  *   (s)  (80) 

1  s      s   L t      t t      t  -l 

A 

*    (s)     =     (K   /J    )  f(J      /J   )s      *     (K   /J   )1     A(s) 

2  pti-sts  tsJ 


(81) 


A         =  as   defined    by   equation     (40) 

B.       CONTROLLER    DESIGN    AND    FREQUENCY    RESPONSE. 

The   problem   of      ship   roll   control    is      a   frequency   domain 
problem:      controlling      against   a      (time      domain)       stochastic 
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disturbance      with      a      generally        known       (frequency      domain) 

spectrum.      Therefore     the  classical   approach   to      control    was 

pursued,    designing     on   the      basis   of      feedback   compensation. 

It      was   assumed      that      roll      angle,    rate,      and      acceleration 

signals   would   be    available   from   shipboard   sensors    (gyros    and 

acceleroraeters) .        The      controller        transfer      function      was 

considered   to    be    of   the    form 

C(s)     =    A(s)/<i>  (s) 

2 
=   K    s      ♦    K    s       +K  (82) 

2  10 

For   purposes   of   block      diagram   development    and   analysis, 

the    following      block    transfer   functions   were      defined,    based 

upon   equations    (79),     (80),    and    (81): 

SS(s)       =     4>    (s)     /  *(s) 

=    K   /J      Ts      ♦     (B   /J    )  s    ♦     (K   /J    )~l/A 

s       s    L  tt  ttJ 

-      (j)         [s      +    2  f  u   s    +    w  1/  A  (83) 

n      *-  t  n  nJ 

ST  (s)       =      *    (s)     /    A  (s) 
2 

=    K   /J       f"(J      /J   )s      ♦     (K    /J   )"]/  A 
p      t    L-     st      s  tsJ 

=    K    /J         |XS       ♦   w    "]/  A  (84) 

p      s      ■-  n-1 

SST  (s)     =      ST(s)     /    SS  (s) 

2  2  2  2 

=    K    /K        fxs      +u    1/Ts      +    2U   s   +   w  (85) 

p      s     L_  n-l    L  t    n  nJ 

G(3)  =         K    /J      Ts       ♦     (B    /J    )  s    ♦     (K    /J    )~1 

s      s   L t      t t      t  rj 

A  -    K    /J       Ps      +     (B    /J    )  s    +     (K   /J    )"] 

s      s    L  t      t  t      t  J 

2         2  2 

w      fs      +    2U   s    +   w  (86) 

a   L. t    n n-1 

2    r  2  ^1 

A   -    uj        s      +    2f  u   s    +  ui 

n    L  t    n  nJ 
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H(s)  =      SST(s)     C(s) 

4 


Where 


=    1/K    TxC    s    +XC    s    ♦  (XC    +wC)s    +    wCs+cjC~| 

sL     2 1  0      n    2 n    1  n   O-l 

2  2 

s      +2fws    +   w  (87) 

t   n  n 


c 

2 

= 

K  K 
P    2 

C 
1 

s 

K  K 
P    1 

c 

0 

3 

K    K 

p  o 

The    "open  loop"   transfer   function,    then,    was 
F(s)  =        G(s)    /      [l    *    G(s)H(s)] 


2 


(K    /J    )  Ts      «■     (B    /J    )  s    ♦    (K    /J    )1 

s      s    L t i t      t  J 

A  (open) 


(88) 


A  (open)    =       A-    (K  /J   )Ts      *    (B   /J   )  s    +    (K  /J   )~|    + 

s      s  *-  tt  ttJ 

-  4  3  2 

♦     (1/J    )  |  XC    s    +  X    C    s    +  (XC    +(K    /J    )C    )s 
s   L     2  1  0        t      t      2 

>(K/J)Cs+(K/J)C~| 

t      t      1  t      t       OJ  (89) 

The   overall,    or   "closed   loop"    response    became 

T(s)  =  *  (s)    /    *   (s) 

F(s)     /      [1    ♦    F(s)] 

r   2  -i 

(K    /J    )|s      ♦    (B    /J    )s    +  (K    /J   )|             (90) 

s      s    L t      t  t      tJ 

A  (closed) 

r  4  3  2 

A(closed)    =    A  +     (1/J    )  I  XC    s    +XC    s    +  (XC    ♦  (K   /J    )C    )s 

+  (K   /J    )C    s    +  (K   /J    )C  (91) 

t      t       1  t      t      O-l 
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The   Characteristic  Equation  of  the   closed  loop   response, 
A(closed) ,  was 


4  3  2 

s       ♦    A    s      ♦    A    s       ♦    A    s      +    A         =0 
3  2  1  0 


J    J 

s    t 

J    J    +J      C 
s   t      st   2 


J   J 

s   t 

J    J    +J      C 
s   t      st   2 


BJ+JB      ♦  J       C 

St         s    t  st    1 

2 

J   J   +J  J   J 

s    t      st  s   t 


(92) 


(93) 


KJ    +    JK    +    BB-2KJ 

st        st        st  tst 

2 
J  J    +J 
s   t      st 


J      C      ♦    K    C 

st    0  t    2 

J   J 

s   t 


(94) 


J    J 

s    t 

J    J    +J      C 
s   t      st   2 


J   J 

s    t 

J    J    +  J      C 

s  t      st   2 


BK+KB      +  KC 

s    t         s    t  t    1 

2 

J   J   -J  J   J 

s   t      st  s   t 


K    K    -    K  +    K    C 
s   t        t  tO 

2 
J   J   -J  J   J 

s   t      st  s   t 

J 


(95) 


(96) 


The  block  diagram  development  is  illustrated  in  Figure 
3-1.  The  first  design  approach  attempted  was  the  modern 
control  theory  state  variable  feedback  method.  The 
difficulty  with  this  method  was  the  inability  to  define  a 
time-domain  objective  or  cost  functional. 

The  next  design  approach  attempted  was  with  the  Computer 
Aided  Design  of  Optimal  Compensators  (CADOC)  program 
[  Anastassakis ,  1977].  This  method  consisted  of  selecting  a 
frequency-domain  cost  functional  and  a  desired  frequency 
response  curve,  and  used  the  function  minimization  procedure 
of  M.J.  Box.  There  were  two  difficulties  encountered  in 
applying  this  program   to   the  problem  of   ship   anti-roll 
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a)  Direct  implementation  of  equations  (79)  and  (82) 


(b)  First  state  of  block  diagram  reduction. 


Figure  3.1. 


Block  Diagram  Development  for  Activated 
Ship/Tank  System 
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(c)  Final  block  diagram. 


(d)  Equivalent  block  diagram  for  servo  theory  analysis 


Figure  3.1. 


Block  Diagram  Development  for  Activated 
Ship/Tank  System 


46 


control.  The  first  problem  was  that  a  portion  of  the 
feedback  path  in  the  physical  system  was  fixed  by 
architecture,  while  the  CADOC  program  required  all 
coefficients  in  the  feedback  path  to  be  adjustable.  The 
second  problem  involved  selecting  an  attainable  frequency 
response  for  CADOC  input.  It  was  concluded  that  the  CADOC 
program  was  best  suited  for  optimizing  a  close  design, 
rather  than  for  initiating  a  design. 

Another  design  approach  attempted  was  in  the  root  locus 
parameter  space.  Two  separate  procedures  exist  here,  and 
both  were  attempted.  The  first  procedure  was  the  standard 
root  Iocjs  approach,  fixing  two  of  the  three  adjustable 
gains  and  varying  the  third.  This  method  amounted  to  an 
exhaustive  search  technique,  attempting  to  locate  the  best 
combination  of  the  three  gains  to  define  system  roots  with 
adequate  lamping. 

The  second  root  locus  parameter  space  procedure  involved 
using  a  parameter  plane  computer  program,  PARAM  A,  which 
exists  at  the  Naval  Postgraduate  School  and  is  available 
from  Dr.  3. J.  Thaler,  but  has  not  been  documented  in  any 
literature  to  his  knowledge.  This  program  varies  two  gains 
which  appear  in  linear  combinations  in  the  coefficients  of 
ths  system  characteristic  equation,  for  a  fixed  value  of  the 
third  gain.  The  output  of  this  program  is  a  two-dimensional 
parameter  plane  showing  families  of  curves  of  constant 
damping  ratio  and  constant  natural  frequency  as  functions  of 
the  two  gains. 

The  problems  associated  with  both  root  locus  parameter 
space  procedures  were  the  exhaustive  search  nature  of  a 
three-variable  system,  and  the  exclusions  of  the  system 
closed-loop  zeroes  from  the  analysis.  This  was  noted  to  be 
a  problem  when  the  gains  found  by  these  procedures  were  used 
in  a  frejuency  response  study  of  the  activated  ship/tank 
system,  and  large  resonant  peaks  were  observed. 
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After  attempting  all  of  the  above  methods  and  abandoning 
them  for  the  reasons  stated,  an  exhaustive  search  of 
closed-loop  frequency  response  curves  was  conducted,  using 
the  program  listed  in  Appendix  C.  An  acceptable  design  was 
accomplished  when  C  =  110,  C  =  -150,  and  C  ■  -80.  This 
corresponded  to  placing  the  compensator  roots  at  s  -  -2.4388 
and  s  =  +0.5638,  which  improved  the  overall  system  damping 
and  provided  the  frequency  response  illustrated  in  Figure 
3-2.  This  response  was  characterized  by  two  low  resonant 
peaks,  one  at  0.582  rad/sec  of  2.85  db  (linear  gain  1.39), 
and  one  at  1.38  rad/sec  of  1.96  db  (linear  gain  1.25).  The 
minimum  response  between  the  two  resonances  occured  at  0.861 
rad/sec,   and    was   -4.75   db    (linear   gain    0.58). 

In  sea  state  6,  the  maximum  waveslopes  at  the  resonant 
frequencies  are  1.44  degrees  at  0.582  rad/sec,  and  8.11 
degrees  at  1.38  rad/sec,  and  the  maximum  waveslope  at  0.861 
rad/sec  is  3.15  degrees.  Using  this  controller  design,  the 
maximum  roll  angles  that  the  ship  should  experience  in  sea 
state  6  would  be  2.0  degrees  at  0.582  rad/sec,  1.83  degrees 
at  0.861  rad/sec,  and  10.1  degrees  at  1.38  rad/sec.  The 
wave  height  power  spectral  density  at  those  frequencies, 
according  to  the  spectra  of  Pierson  and  Moskowitz,  are  205 
ft2  /Hz  at  0.582  rad/sec  (0.093  Hz),  80  ft 2  /Hz  at  0.861 
rad/sec  (0.137  Hz),  and  12  ft2/Hz  at  1.38  rad/sec  (0.219 
Hz) .  Therefore  it  was  concluded  that  even  though  the  pure 
numbers  indicated  a  poor  roll  control  at  the  higher 
frequency  resonant  peak,  the  waves  at  that  frequency  lack 
sufficient  power  to  cause  the  ship  to  attain  the  maximum 
theoretical   roll   angle. 
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FREQUENCY  RESPONSE,  ACTIVATED  TANKS 
oCLQSED  LOOP  RESPONSE,  MRG  (DB)  VS.  LOG  N 


-0.80 


0.80 


Figure  3-2. 


Frequency  Response  of  Activated  Ship/Tank 
System,  Proportional  Plus  Rate  Plus  Acceleration 
Control 
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IV.   EFFECTS  OF  VESSEL  DISPLACEMENT  VARIATIONS 

The  analysis  of  ship  motions  and  roll  control  thus  far 
had  been  conducted  under  the  assumption  that  the 
displacement  of  the  vessel  was  constant.  Over  the  short 
term,  the  assumption  may  be  true  enough,  but  for  a  ship 
underway,  consuming  fuel  and  storas,  the  assumption  cannot 
hold  true  beyond  a  few  days  at  best. 

As  a  ship's  displacement  decreases,  the  underwater 
surface  area  of  the  hull  decreases  also,  in  keeping  with  the 
Principle  of  Archimedes.  This  in  turn  leads  to  a  lowering 
of  the  Canter  of  Buoyancy,  which  results  in  a  lowering  of 
the  Metacenter.  The  effect  on  the  metacentric  height, 
however,  is  not  so  clear,  as  the  exact  location  of  the 
Center  of  Gravity  in  the  ship  depends  more  on  the  placement 
of  the  onboard  weight  than  on  the  total  displacement.  In 
general,  fuel  tanks  are  located  low  in  a  ship;  thus 
consumption  of  fuel  is  removal  of  low  weight,  which  shifts 
the  Center  of  Gravity  upward,  decreasing  the  metacentric 
height.  It  should  be  noted  here  that  ballasting  empty  fuel 
tanks  with  sea  water  has  the  effect  of  restoring  the  ship  as 
muoh  as  possible  to  the  original  condition  of  statical 
stability  which  existed  at  the  beginning  of  the  voyage,  and 
extends  the  time  frame  of  applicability  of  the  assumption  of 
constant  displacement. 

Varying  the  ship's  displacement  therefore  varies  the 
natural  frequency  of  the  ship,  and  effectively  de-tunes  the 
anti-roll  tank  system.  The  only  physical  adjustment  that 
can  be  made  to  a  passive  tank  installation  underway  is  to 
vary  the  quantity  of  water  in  the  U-tube.  This  has  the 
effect  of  changing  the  limits  of  integration  of  the  design 
equations,  and  therefore  changes  the  parameters  of  the  tank. 
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Additional  adjustments  are  available  on  the  activated 
tank  system.  These  are  the  gain  of  the  pump,  which  is  a 
function  of  the  pump  rotational  speed,  and  the  feedback 
gains  from  the  sensors  to  the  pump  controls.  Thus,  if  a 
method  exists  for  determining  on-line  the  various  ship 
parameters  as  they  change,  the  controller  for  the  activated 
tank  system  can  be  made  adaptive.  Various  schemes  may  be 
considered,  from  Kalman  filter  type  parameter  estimators  to 
FFr  devices,  and  this  field  is  left  as  a  topic  for  further 
possible  research. 
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V.   CONCLUSIONS 

In  studying  ship  roll  control,  it  was  necessary  to 
become  familiar  with  the  nature  of  ocean  waves.  The 
waveslope  spectrum  of  the  ocean  which  was  generated,  when 
used  in  conjunction  with  the  Pierson-Moskowitz  wave  power 
spectra,  gave  much  insight  into  the  problem.  It  was  learned 
that  the  sea  surface  does  not  and  can  not  assume  the  slopes 
that  an  observer  on  a  ship  often  estimates. 

A  typical  ship  on  the  ocean  magnifies  the  waveslopes 
near  the  resonant  frequency  by  a  factor  of  about  15.  By 
installing  passive  anti-roll  tanks,  this  waveslope 
magnification  factor  was  split  into  two  resonant  peaks,  one 
with  a  factor  of  5,  and  the  other  with  a  factor  of  3. 
Activating  the  tank  system  further  decreased  the  waveslope 
magnification  factors  to  1.4  and  1.3,  respectively. 

Several  assumptions  have  been  made  over  the  years  in  the 
study  of  roll  control  by  U-tube  tank  systems,  and  some  of 
them  were  sources  of  problems  in  the  design  of  physical 
systems.  The  first  assumption  was  that  a  ship  has  constant 
parameters:  displacement,  metacentric  height,  and  natural 
period.  This  assumption  holds  true  only  if  the  ship  does  not 
consume  any  fuel  or  stores,  or  if  these  weights  are  replaced 
with  appropriate  quantities  of  suitably  placed  ballast. 

Another  assumption  was  that  the  tanks  could  be  placed  in 
the  ship  in  such  fashion  as  to  drive  the  mutual -coupling 
inertia  coefficient  to  zero.  This  assumption  had  two 
limitations.  The  first  pitfall  was  the  direct  connection 
between  this  assumption  and  the  assumption  of  constant  ship 
parameters.  The  second  pitfall  was  that  the  tanks  had  to  be 
placed  in  an  impossible  location  on  the  ship  to  accomplish 
the  objective  of  the  assumption. 
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Roll  control  by  pump  activated  tanks  is  perhaps  the  best 
form  of  roll  control  for  slow  moving  vessels,  or  for  vessels 
conducting  operations  either  at  slow  speed  or  stopped,  where 
activated  fin  stabilizers  lose  their  effectiveness.  It 
should  be  noted  that  a  pump  activated  tank  system  can  be 
used  to  cause  rolling,  as  an  aid  for  breaking  out  of  ice. 
It  can  also  be  used  to  correct  list,  including  damage-caused 
list,  as  a  safer  alternative  to  counter- flooding.  Therefore 
the  pump  activated  U-tube  tank  system  is  the  most  versatile 
form  of  coll  control  that  a  multi-purpose  vessel  can  have 
installed. 
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APPENDIX  A 


Program  Listing  for  Frequency  Response  of 
Unmodified  Ship,  with 
FORTRAN  Subroutine  PARAM 


The  Fcsquency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IBM 
Systera/360  machine  at  the  W.  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
manual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  8 


Program  Listing  for  Frequency  Response  of 

Ship  with  Passive  Tanks,  with 

FORTRAN  Subroutine  CALC 


The  Frequency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IBM 
System/360  machine  at  the  W.  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
manual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  C 


Program  Listing  for  Frequency  Response  of 
Activated  Ship/TanJc  System,  with 
FORTRAN  Subroutine  CALCA 


The  Frequency  Response  program  was  written  in  Digital 
Simulation  Language  (DSL/360)  to  be  run  on  the  IBM 
System/360  machine  at  the  W.  R.  Church  Computer  Center, 
Naval  Postgraduate  School.  To  use  the  program,  a  deck  of 
plot  identification  cards  must  be  appended  following  the 
FORTRAN  Subroutine,  in  accordance  with  the  DSL/360  user's 
manual. 

The  labels  and  variable  names  used  in  the  program  are 
defined  in  the  Glossary  for  Frequency  Response  Programs, 
Appendix  D  of  this  thesis. 
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APPENDIX  D 


Slossary  for  Frequency  Response  Programs 


MAIN 


SYMBOL 

AO, A1 , A2, A3 

BS 

BT 

CLOOP 

COMP 

D 

G 

GIMAG 

GREAL 

H 

IM 

IMH 

JS 

JSI 

JT 

KS 

KT 

LOGW 

HAG 

MA3H 

MDB 

MDBH 

N 

OLOOP 

PHI 

MEANING 

Coefficients  of  characteristic  equation 

Damping  coefficient  of  ship 

Overall  damping  coefficient  of  tank-fluid 

Closed-loop    transfer    function 

Compensator    (controller)    transfer   function 

Denominator      polynomial   of     open-loop    transfer 

function 

Open-loop  transfer  function 

Imaginary  part  of  open-loop  transfer  function 

Real  part  of  open-loop  transfer  function 

Closed-loop  transfer  function 

Imaginary  part  of  open-loop  transfer  function 

Imaginary  part  of  open-loop  transfer  function 

Roll  inertia  coefficient  of  ship 

Mutual-coupling  inertia  coefficient 

Equivalent  inertia  coefficient  of  tank-fluid 

Static  righting  coefficient  of  ship 

Equivalent  static  coefficient  of  tank-fluid 

Natural  logarithm  of  radian  frequency 

Magnitude  of  open-loop  transfer  function 

Magnitude  of  closed-loop  transfer  function 

Decibel  gain  of  open-loop  transfer  function 

Decibel  gain  of  closed-loop  transfer  function 

Numerator  polynomial  of  transfer  function 

Open-loop  transfer  function 

Phase  of  open-loop  transfer  function 
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Q  Denominator   polynomial  of  open-loop  transfer 

function 
S  LaPlace  transform  variable 

SST  Ship-tank  response  in  feedback  path 


SUBROUTINES 

SYMBOL         MEANING 

A  Cross-sectional  area  of  vertical  leg  of  U-tube 

AS  Cross-sactional  ar=a  of   horizontal   leg   of 

U-tube 
BEAM  Transverse  dimension  of  ship  (beam) 

BP  Velocity  head-loss  coefficient  of  pump 

BS  Damping  coefficient  of   ship 

BT  Overall  damping  coefficient  of  tank-fluid 

BT1  Skin-friction   damping   coefficient   of   tank 

fluid 
D  Perpendicular  distance  from  center  of  rotation 

of  ship  to  horizontal  cross-leg  of  U-tube 
DABS  Absolute  value  (magnitude)  of  DCAL 

DCAL  Calculated   value    of    D   to   force    the 

mutual-coupling  inertia  coefficient  to  zero 
DISP  Displacement  of  ship 

E1,E2  Incremental   diffenance  between   current   and 

previously   calculated    values   of    ship 

parameters 
GM  Metacentric  height 

GRAV  Acceleration  of  gravity 

H  Height  of   fluid  column   in   vertical  leg   of 

U-tube 
11,12,13,14    Evaluation   of  integrals  as  described   in  main 

body  of  thesis 
JS  Roll  inertia  coefficient  of  ship 
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JS1  Roll  inertia  coefficient  of  unmodified  ship 

JST  Mutual-coupling  inertia  coefficient 

JT  Equivalent  inertia  coefficient  of  tank-fluid 

KRAD  Radius  of  gyration  of  the  mass  of  the  ship 

KS  Static  righting  coefficient  of  ship 

KS1  Static  righting  coefficient  of  unmodified  ship 

KT  Equivalent   static   coefficient   of   tank-fluid 

L  Longitudinal      dimension   of     a   vertical      leg   of 

the   U-tube 

LEN  Longitudinal    dimension   of   ship    (length) 

LMDA  Design   ratio   X 

N  Number   of   iterations   of    design    procedure 

R  Lever   arm   of    U-tube 

RHO  Mass   density    of   tank-fluid 

RHDG  Weight   density    of   tank-fluid 

TJ  Temporary   JS 

TK  Temporary    KS 

TPHI  Natural   period   of  roll   of   ship 

WIDE  Transverse   dimension   of    vertical    leg   of    U-tube 

WN2  Natural   frequency   of   ship-tank   system,    squared 

X  Design   ratio    X 
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APPENDIX    E 
Waveslope   Spectrum   of    the   Ocean 

The  wavaslope  spectrum  of  the  ocean  was  generated  from 
the  relationships  of  the  trochoidal  wave  model.  It  is 
intended        to  be        used  in        conjunction  with        the 

Pierson-Moskowitz   wave    power   spectra. 

This  appendix  includes  the  listing  of  the  DSL/360 
prDgram  which  generated  the  spectrum,  and  listings  of  the 
spectra   foe    sea   states    3,    4,    5,    and   6. 
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WAVESLOPE    SPECTRUM    OF    THE    OCEAN 
SIGNIFICANT    HEIGHT    =    3.3    FEET 
******** (SEA    STATE    3)  ******** 

W  PHIMAX 

1.OO00E-01  9.3579E-03 

1.0593E-01  1.0500E-02 

1.  1220E-01  1.  1781E-02 

1.1885E-01  1.3218E-02 

1.2589E-01  1.4831E-02 

1.3335E-01  1.6641E-02 

1.4125E-01  1.8671E-02 

1.4962E-01  2.0950E-02 

1.5849E-01  2.35C6E-02 

1.6788E-01  2.6374E-02 

1.7783E-01  2.9592E-02 

1.8836E-01  3.3203E-02 

1.9953E-01  3.7254E-02 

2.1135E-01  4.1800E-02 

2. 2387E-01  4.6900E-02 

2.3714E-01  5.2623E-02 

2.5119E-01  5.9044E-02 

2.6607E-01  6.6249E-02 

2.8184E-01  7.4332E-02 

2.9854E-01  8.3402E-02 

3.  1623E-01  9.3579E-02 

3.3497E-01  1.0500E-01 

3.5481E-01  1.1781E-01 

3.7584E-01  1.3218E-01 

3.9811E-01  1.4831E-01 

4.2170E-01  1.6641E-01 

4.4668E-01  1.8671E-01 

4.7315E-01  2.0950E-01 

5.0119E-01  2.3506E-01 

5.3088E-01  2.6374E-01 

5.6234E-01  2.9592E-01 

5.9566E-01  3.3203E-01 

6.3096E-01  3.7254E-01 

6.6834E-01  4.1800E-01 

7.0794E-01  4.6900E-01 

7.4989E-01  5.2623E-01 

7.9433E-01  5.9044E-01 

8.4139E-01  6.6248E-01 

8.9125E-01  7.4332E-01 

9.4406E-01  8.3402E-01 

1.0000E    00  9.3578E-01 

1.0593E    00  1.0500E    00 

1.  1220E    00  1. 1781E    00 

1.1885E    00  1.3218E    00 

1.2589E    00  1.4831E    00 

1.3335E    00  1.6641E    00 

1.4125E    00  1.8671E    00 

1.4962E    00  2.0949E    00 

1.5849E    00  2.35C6E    00 

1.6788E    00  2.637UE    00 

1.7783E    00  2.9592E    00 

1.8836E    00  3.3203E    00 

1.9953E    00  3.7254E    00 
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WAVESLOPE  SPECTRUM  OF  THE  OCEAN 

SIGNIFICANT  HEIGHT  =    3.3  FEET 
******** (SE A  STATE  3)  ******** 

tf  PHIMAX 

2.  1135E  00  4.1800E  00 
2.2387E  00  4.6900E  00 
2.3714E  00  5.2623E  00 
2.5119E  00  5.9044E  00 
2.6607E  00  6.6248E  00 
2.8184E  00  7.4331E  00 
2.9854E    00  8.1851E  00 

3.  1623E  00  8. 1851E  00 
3.3496E  00  8.1851E  00 
3.5481E  00  8. 1851E  00 
3.7584E  00  8.1851E  00 
3.9811E  00  8. 1851E  00 
4.2169E  00  8.1851E  00 
4.4668E  00  8. 1851E  00 
4.7315E  00  8.1851E  00 
5.0118E  00  8.1851E  00 
5.3088E  00  8. 1851E  00 
5.6234E  00  8.1851E  CO 
5.9566E  00  8. 1851E  00 
6.3095E    00  8.1851E  00 


74 


WAVESL3PE    SPECTRUM    OF    THE    OCEAN 
SIGNIFICANT    HEIGHT    =6.9    FEET 
********(SEA    STATE    4)  ******** 

W  PHIMAX 

1.0000E-01  1.9566E-02 

1.0593E-01  2.1954E-02 

1.1220E-01  2.4633E-02 

1.1885E-01  2.7638E-02 

1.2589E-01  3.1011E-02 

1.3335E-01  3.4795E-02 

1.4125E-01  3.9040E-02 

1.4962E-01  4.3804E-02 

1.5849E-01  4.9149E-02 

1.6788E-01  5.5146E-02 

1.7783E-01  6.1875E-02 

1.8836E-01  6.9424E-02 

1.9953E-01  7.78S5E-02 

2.1135E-01  8.7400E-02 

2.2387E-01  9.8065E-02 

2.3714E-01  1.1003E-01 

2.5119E-01  1.2346E-01 

2.6607E-01  1.3852E-01 

2.8184E-01  1.5542E-01 

2.9854E-01  1.7439E-01 

3.  1623E-01  1.9566E-01 

3.3497E-01  2.1954E-01 

3.5481E-01  2.4633E-01 

3.7584E-01  2.7638E-01 

3.9811E-01  3.1011E-01 

4.2170E-01  3.4795E-01 

4.4668E-01  3.9040E-01 

4.7315E-01  4.3804E-01 

5.0119E-01  4.9149E-01 

5.3088E-01  5.5146E-01 

5.6234E-01  6.1874E-C1 

5.9566E-01  6.9424E-01 

6.3096E-01  7.78S5E-01 

6.6834E-01  8.7400E-01 

7.0794E-01  9.8064E-01 

7.4989E-01  1. 1003E    00 

7.9433E-01  1.2346S    00 

8.4139E-01  1.3852E    00 

8.9125E-01  1.5542E    00 

9.4406E-01  1.7439E    00 

1.0000E    00  1.9566E    00 

1.0593E    00  2.1954E    00 

1.  1220E    00  2.4633E    CO 

1.1885E    00  2.7638E    00 

1.2589E    00  3.1011E    00 

1.3335E    00  3.4794E    00 

1.4125E    00  3.9040E    00 

1.4962E    00  4.3803E    00 

1.5849E    00  4.9148E    00 

1.6788E    00  5.5145E    00 

1.7783E    00  6.  1874E    00 

1.8836E    00  6.9424E    00 

1.9953E    00  7.78S5E    00 
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WAVESLOPE  SPECTRUM  OF  THE  OCEAN 
SIGNIFICANT  HEIGHT  =6.9  FEET 
****#***(SEA  STATE  4) ******** 

W  PHIMAX 

2.  1135E  00  8.1851E  00 
2.2387E  00  8.1851E  00 
2.3714E  00  8.1851E  00 
2. 5119E  00  8. 1851E  00 
2.6607E  00  8.1851E  00 
2.8184E  00  8.1851E  00 
2.9854E    00  8.1851E    00 

3.  1623E  00  8.  18512  00 
3.3496E  00  8. 1851E  00 
3.5481E  00  8.1851E  CO 
3.7584E  00  8.1851E  00 
3.9811E  00  8. 1851E  00 
4.2169E  00  8.1851E  00 
4.4668E  00  8.1851E  00 
4.7315E  00  8.1851E  00 
5.01 18E  00  8. 1851E  00 
5.3088E  00  8.1851E  00 
5.6234E  00  8.1851E  00 
5.9566E  00  8. 1851E  00 
6.3095E    00  8.1851E    CO 
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tfAVESLOPE    SPECTBUM    OF    THE    OCEAN 
SIGNIFICANT    HEIGHT    =    10.     FEET 
********(SEA    STATE    5)  ******** 

H  PHIMAX 

1.0000E-01  2.8357E-02 

1.0593E-01  3.1817E-02 

1.1220E-01  3.5700E-02 

1.1885E-01  4.0056E-02 

1.2589E-0  1  4.4943E-02 

1.3335E-01  5.0427E-02 

1.4125E-01  5.6580E-02 

1.4962E-01  6.3484E-02 

1.5849E-01  7.1230E-02 

1.6788E-01  7.9921E-02 

1.7783E-01  8.9673E-02 

1.8836E-01  1.0062E-01 

1.9953E-01  1.1289E-01 

2.1135E-01  1.2667E-01 

2.2387E-01  1.4212E-01 

2.3714E-01  1.5946E-01 

2.5119E-01  1.7892E-01 

2.6607E-01  2.0075E-01 

2.8184E-01  2.2525E-01 

2.9854E-01  2.5273E-01 

3.1623E-01  2.8357E-01 

3.3497E-01  3.1817E-01 

3.5481E-01  3.5699E-01 

3.7584E-01  4.0056E-01 

3.9811E-01  4.4943E-01 

4.2170E-01  5.0427E-01 

4.4668E-01  5.6580E-01 

4.7315E-01  6.3484E-01 

5.0119E-01  7.1230E-01 


5.3088E-01  7.9921E-01 

5.6234E-01  8.9673E-01 

5.9566E-01  1.0061E    00 

6.3096E-01  1.1289E    00 

6.6834E-01  1.2667E    00 

7.0794E-01  1.4212E    00 

7.4989E-01  1.5946E    00 

7.9433E-01  1.7892E    00 

8.4139E-01  2.0075E    00 

!    8.9125E-01  2.2525E    00 

9.4406E-01  2.5273E    00 

1.0000E    00  2.8357E    00 

1.0593E    00  3.1817E    00 

1.  1220E    00  3.5699E    00 

1.1885E    00  4.0055E    00 

1.2589E    00  4.4943E    00 

1.3335E    00  5.0427E    00 

1.4125E    00  5.6580E    00 

1.4962E    00  6.3483E    00 

1.5849E    00  7.1230E    00 

1.6788E    00  7.9921E    00 

1.7783E    00  8. 1851E    00 

1.8836E    00  8.1851E    00 

1.9953E    00  8.1851E    CO 
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WAVESLOPE    SPECTRUM   OF    THE    OCEAN 

SIGNIFICANT    1 

iEIGHT    =    10. 

FEET 

******** (SEA 

STATE   5)  ***3 
PHI! 

***** 

W 

1AX 

2.  1135E 

00 

8.1851E 

00 

2.2387E 

00 

8.1851E 

00 

2.3714E 

00 

8. 1851E 

00 

2.5119E 

00 

8.  1851E 

00 

2.66  07E 

00 

8.1851E 

00 

2.8184E 

00 

8.  1851E 

CO 

2.9854E 

00 

8.1851E 

00 

3.  1623E 

00 

8.  1851E 

00 

3.3496E 

00 

8.1851E 

00 

3.5481E 

00 

8.  1851E 

00 

3.7584E 

00 

8.  1851E 

00 

3.9811E 

00 

8.  1851E 

00 

4.2169E 

00 

8.  1851E 

00 

4.4668E 

00 

8.  1851E 

00 

4.7315E 

00 

8.1851E 

00 

5.0118E 

00 

8.  1851E 

00 

5.3088E 

00 

8.1851E 

00 

5.6234E 

00 

8.1851E 

CO 

5.9566E 

00 

8.  1851E 

00 

6.3095E 

00 

8.1851E 

00 
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WAVESLDPE    SPECTEUM    OF    THE    OCEAN 
SIGNIFICANT    HEIGHT    =    15.     FEET 
********(SEA    STATE    6)  ******** 

W  FHIMAX 

1.0000E-01  4.2536E-02 

1.0593E-01  4.7726E-02 

1.1220E-01  5.3549E-02 

1.1885E-01  6.0083E-02 

1.2589E-01  6.7415E-02 

1.3335E-01  7.5641E-02 

1.4125E-01  8.4870E-02 

1.4962E-01  9.5226E-02 

1.5849E-01  1.0685E-01 

1.6788E-01  1.1988E-01 

1.7783E-01  1.3451E-01 

1.8836E-01  1.5092E-01 

1.9953E-01  1.6934E-01 

2.1135E-01  1.9000E-01 

2.2387E-01  2.1318E-01 

2.3714E-01  2.3920E-01 

2.5119E-01  2.6838E-01 

2.6607E-01  3.0113E-01 

2.8184E-01  3.3787E-01 

2.9854E-01  3.7910E-01 

3.  1623E-01  4.2536E-01 

3.3497E-01  4.7726E-01 

3.5481E-01  5.3549E-01 

3.7584E-01  6.0083E-01 

3.9811E-01  6.7414E-01 

,   4.2170E-01  7.5640E-01 

4.4668E-01  8.4870E-01 

4.7315E-01  9.5225E-01 

5.0119E-01  1.0684E    CO 

5.3088E-01  1.1988E    00 

5.6234E-01  1.3451E    00 

■    5.9566E-01  1.5092E    00 

6.  3096E-01  1.6934E    00 

6.6834E-01  1.9000E    00 

7.0794E-01  2.1318E    00 

7.4989E-01  2.3920E    00 

7.9433E-01  2.6838E    00 

8.4139E-01  3.0113E    00 

8.9125E-01  3.3787E    CO 

9.4406E-01  3.7910E    00 

1.0000E    00  4.2536E    00 

1.0593E    00  4.7726E    00 

1.  1220E    00  5.3549E    00 

1.1885E    00  6.0083E    00 

1.2589E    00  6.7414E    00 

1.3335E    00  7.5640E    00 

1.4125E    00  8.1851E    00 

1.4962E    00  8.1851E    00 

1.5849E    00  8.1851E    00 

1  .6788E    00  8.1851E    00 

1.7783E    00  8. 1851E    00 

1.8836E    00  8. 1851E    00 

1.9953E    00  8. 1851E    00 
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WAVESLOPE    SPECTRUM    OF    TEE    OCEAN 

SIGNIFICANT    HEIGHT    =    15.  FEET 
******** ( SEA    STATE    6) ******** 

8  PHIMAX 

2.1135E    00  8.1851E  00 

2.2387E    00  8.1851E  CO 

2.3714E    00  8.1851E  00 

2.5119E    00  8. 1851E  00 

2.6607E    00  8.1851E  00 

2.8184E    00  8.1851E  00 

2.9854E    00  8. 1851E  00 

3.  1623E    00  8. 1851E  00 

3.3496E    00  8. 1851E  00 

3.5481E    00  8. 1851E  CO 

3.7584E    00  8.1851E  00 

3.9811E    00  8.1851E  00 

4.2169E    00  8.1851E  00 

4.4668E    00  8. 1851E  00 

4.7315E    00  8.1851E  00 

5.0118E    00  8. 1851E  00 

5.3088E    00  8.1851E  00 

5.6234E    00  8.1851E  00 

5.9566E    00  8. 1851E  00 

6.3095E    00  8.1851E  CO 
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